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[57] ABSTRACT 

An apparatus and computer-implemented polygon bounding 
method generate, in a first step, a binary array representation 
of a closed contour input object relative to a reference 
coordinate system. Subsequently, in a single pass through 
the binary array, a bounding polygon is generated as a 
sequence of parallel "bands," each representing a pair of 
opposing sides of the polygon. The polygon bands are 
oriented at successive generally uniform rotations of the 
reference coordinate system and are positioned to define the 
bounds of the closed contour for each rotation angle. The 
spacing and positioning of the bands can thus be used to 
define polygon size and position at each angle. In this way, 
comparisons between a plurality of polygons can be made 
systematically using defined geometrical criteria. 

18 Claims, 4 Drawing Sheets 
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APPARATUS AND FAST METHOD FOR 
CONVEX-POLYGON BOUNDING OF 
CLOSED CONTOURS 

This application is a continuation of application Ser. No. 
08/425^14, filed Apr. 17, 1995, now abandoned, which is a 
continuation of application Ser. No. 08/036,663, filed Mar. 
24, 1993, now abandoned. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to closed contour 
object characterization, and more particularly to the isola- 
tion and classification of closed contour objects using 
convex-polygon bounding structures. 

It is sometimes convenient to isolate and characterize 
objects in two dimensional and three dimensional space by 
enclosing them within a defined polygon (or polyhedron in 
3D) whose size and shape provide information about the 
objects they enclose. For example, bounding polygons can 
be used to distinguish the components of an array of objects 
and determine rough approximations of their shape. Specific 
applications in which polygon bounding has proven useful 
include pattern recognition, in the case of two dimensional 
objects, and virtual imaging in the case of objects defined in 
three dimensional space. 

In the field of pattern recognition, a convex-polygon 
bounding method has been proposed by K. Yamamoto et al, 
Recognition of Hand-Printed Characters By An Outermost 
Method, Pattern Recognition, Vol. 12, pp. 229-236 (1980). 
The Yamamoto et al method approximates the convex hull 
of a discrete set of points. Approximate convex hull struc- 
tures are generated by identifying, in the object to be 
bounded, a series of "outermost" contour points. Those 
points are sequentially connected to form the sides of the 
bounding polygon. For example, the letter "A" might be 
bounded by a polygon having a short horizontal upper side, 
a long horizontal lower side, and a pair of slanted interme- 
diate sides. 

Each side of the polygon so generated is inherently 
oriented at some angular rotation with respect to a fixed 
frame of reference. The angular rotations are determined 
solely by the positions of the contour points, which in turn 
are determined by the shape of the closed contour to be 
bounded. Thus, in the prior method, the angular rotations of 
the polygon sides, and even the number of sides, typically 
differ from polygon to polygon when different contours are 
considered. For example, the sides of a polygon constructed 
to bound the letter "A" would generally be fewer in number 
and have different angular rotations than a polygon bound- 
ing the letter "B". Polygons generated for different contours 
in accordance with the Yamamoto et al method therefore 
have few, if any, characteristics in common, and cannot be 
easily compared. 

It is submitted that contour analysis could be significantly 
improved using convex-polygons having common 
characteristics, such as a predetermined number of sides and 
uniform rotational characteristics. In that way, information 
such as polygon size with respect to selected angular rota- 
tions could be easily compared relative to a plurality of 
polygons, regardless of the shapes of the bounded contours. 
It would be further desirable to provide a method for 
generating such convex-polygons in an efficient manner 
using minimal processing resources. For example, in known 
character recognition systems, closed contours are detected 
during one pass over a binary image array representation. It 
would be convenient to obtain at the same time a convex- 
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polygon bound on that contour. This would simplify clas- 
sification of the object represented, as well as the compu- 
tation of intersections, point membership and other 
topological relations. 

5 What is needed, therefore, is a class of bounding polygons 
allowing computational advantages that are inherent in their 
geometrical characteristics, as well as a fast method for 
generating such polygons as an extension to the classical 
one-pass algorithm for finding connected components in 

30 binary images. 

SUMMARY OF THE INVENTION 

In accordance with the foregoing objectives, an apparatus 

1S and computer-implemented polygon bounding method 
generate, in a first step, a binary array representation of a 
closed contour input object relative to a reference coordinate 
system. Subsequently, in a single pass through the binary 
array, a bounding polygon is generated as a sequence of 

2Q parallel "bands," each representing a pair of opposing sides 
of the polygon. The polygon bands are oriented at successive 
generally uniform rotations of the reference coordinate 
system and are positioned to define the bounds of the closed 
contour for each rotation angle. The spacing and positioning 

^ of the bands can thus be used to define polygon size and 
position at each angle. In this way, comparisons between a 
plurality of polygons can be made systematically using 
defined geometrical criteria. 

In a preferred aspect of the invention, each polygon band 

30 pair is defined by a pair of contour bounding displacements. 
The contour bounding displacements represent minimum 
and maximum contour distance values with respect to a 
selected axis (e.g. the horizontal "x" axis) of the reference 
coordinate system, as it is positioned through each angle of 

35 rotation. Advantageously, the contour bounding displace- 
ments at each rotation angle can be determined by mapping 
the contour displacement values into the coordinate system 
at that angle. Contour point displacements (interval end 
points) are determined for each row of the input array, for 

40 each rotation angle. Global minimum and maximum dis- 
placements for each rotation angle are found by replacing 
the largest value found at a given angle by the next value 
computed, if that next value is still larger. The final mini- 
mum and maximum for each angle are stored in an output 

as table. 

In still another aspect of the invention, the mapping 
procedure is performed by selecting rotation angle approxi- 
mation values representing rational tangents of the form 
me/no, wherein successive m e and successive n 6 values 
50 differ by at most 1. This facilitates determination of the 
contour displacement, for each, rotation angle, by simply 
incrementing and decrementing the reference system dis- 
placements in accordance with a predetermined sequence 
over the range of angles. 

55 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram showing a data processing 
system constructed in accordance with the present invention. 
60 FIG. 2 is a flow diagram showing a preferred method in 
accordance with the present invention. 

FIG. 3a is a diagrammatic illustration of a document 
image containing a closed contour object. 

FIG. 3b is a diagrammatic illustration of a binary array 
65 representation of the document image showing a reference 
coordinate system relative thereto together with selected 
incremental rotations of the reference coordinate system 
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(shown in phantom), and further showing minimum and [x^y), x^y)], where y ranges between the first scan line of 

maximum contour bounding displacements relative to the contour, y h (bottom), and the line at which the contour 

selected rotations of the reference coordinate system. closes, y, (top). 

FIG. 4 is a diagrammatic illustration of an interval table If only a rudimentary polygon is desired, such as a 

generated in accordance with the present invention. s bounding rectangle, the reference coordinate system scan 

FIG. 5 is a flow diagram showing a fast computational mtc ™\ s ^ ucncc ^ * <* d " cc % ^ * bounding 

method for performing contour point interval mapping. ™ Un &> fo / c ™"P k ' le * Slde 15 ™*{*iW> the 

„ ■ - . • s *dc is maxlxJy)}, and the bottom and top are y b and y„ 

FIG. 6 is a detailed diagrammatxc illustration of a binary ^^ely. A disadvantage of bounding rectangles and 

array representation of a closed contour image showing a JQ oAcj ^ le sfa fc ^ ^ m rf ^ {q[ ^ 

relationship between minimum and maximum contour dis- classification> ^ for determining the actual proximity of 

placement values in successive scan lines warranting a lura , dosed rantour objec(s ^ ^ becanse tbe boujjdi 

termination of further rotational processing for the first of rectangl( , ^ genera]ly a pQor approxinlatioil of Ae 

the illustrated scan line pair. in order to improve object classification and proximity 

DETAILED DESCRIPTION OF THE 15 determination, tne notion of a bounding rectangle must be 

PREFERRED EMBODIMENTS generalized to that of a bounding convex polygon consisting 

of an arbitrary fixed number of pairs of opposite parallel 

Turning now to FIG. 1, a character bounding system for sides, instead of just two such pairs. This requires that the 

evaluating closed contour images in accordance with a bounds of the object be computed along multiple fixed 

preferred embodiment of the invention is constructed from 20 directions in a plane. 

a data processing apparatus 10. The data processing appa- As shown in FIG. 3fc, it is desired to construct a convex 

ratus 10 includes a processing unit 12, a data input system polygon having opposing sides or bands at generally uni- 

14 for the input of data as well as command and control form incremental rotations of the reference coordinate sys- 

information, a data output system 16, and one or more data tern 36. This is done by mapping the scanline interval end 

storage resources 18. "The processing unit 12 includes a 25 displacement values from the reference coordinate system 

plurality of storage register locations 20 as well as local through a sequence of angular rotations. Because the bounds 

program memory 22 for storing and processing instructions on a contour at an angle 6 are also the bounds at 8+ji, the 

from a compiled computer program 24. The data processing rotational domain can be restricted to O^0<n. It is assumed 

apparatus 10 may be constructed from any number of that x^yJQ) and x/y,6) are the minimum/maximum limits of 

programmable general purpose computers and data process- 30 a scan interval of the reference coordinate system 36, rotated 

ing systems. Preferably, the apparatus includes a document by angle 6. The coordinate system 36 is assumed to have its 

scanner as part of the data input system 14, or some other origin at the lower left hand corner of the image. To maintain 

mechanism for inputting data files representing input consistency, y b becomes x r and y r becomes x 1 under a 

images. The data processing apparatus 10 is configured to rotation left by ^/2, with the mapping changed accordingly 

perform the method of the present invention using a soft- 35 for a right rotation. 

ware program written in the well-known "C" programming Referring now to FIG. 2, the first step in the polygon 

language, although it will be understood that many other construction process following input of the array 34 is to 

programming languages could be used to generate a suitable initialize in process step 38 an interval table 40, as shown in 

instruction set. FIG. 4. The interval table 40 is a data output structure that 

CONSTRUCTION 40 defines for each angular rotation of the reference coordinate 

Referring now to FIG. 2, a method for generating closed system 36, the global minimum and maximum contour 

contour bounding polygons is shown diagrammatically. In a bounding intervals x J (y,G), x r (y,0), and the scanline y on 

first process step 30, an input graphics image 32, as shown which they are found. The interval table 40 thus provides the 

in FIG. 3a, having a closed contour, is converted using a output representation of the constructed bounding polygon, 

document image scanner of the data input system 14, to a 45 Following initialization of the interval table 40, a first line of 

binary array representation 34, as shown in FIG. 36. The the input array 34 is scanned in process step 42 and a contour 

binary array representation 34 is a data structure defining a bounding interval x^y^B^ and xj(y nt d^) is determined for 

plurality of array locations arranged in a series of rows and the reference coordinate system 36 (0„=O) in process step 44. 

columns relative to a reference coordinate system 36. By In process step 46, the left and right bounding interval 

convention, array locations having a value of 1 are assumed so displacements x 1 (y n ,0 w ) and x^y^G,,) are separately com- 

to identify darkened portions of the closed contour image, pared against the corresponding 6„ entries for the previous 

while array locations having 0*s represent undarkened por- scanline (if any) in the interval table 40 (i.e.,Xi(y B -i,B), 

tions of the scanned image. As an alternative to scanning, a X-Ay„-ifi)- If one or both are larger, they are stored as a new 

previously generated array representation 34 could be input local minimum, maximum bounding interval 9„. In process 

to the data processing apparatus 10 as a data file by con- 55 step 48, the coordinate system rotation is implemented by 

ventional file handling methods. mapping the interval values x 1 (y,0 l1 ), x£y,Q^ to the next 

The preferred method for generating a bounding polygon rotation angle 6 n+1 . In step 50, the process returns to step 44 

representation of the closed contour input image 32, is to to determine the next left and right interval displacements 

scan the binary array representation 34 in a line -by-line Xi(y>9m.i)» x r (y,0 n+1 ). Process steps 44-50 are repeated until 

(row-by-row) sequence in the direction of the "y" axis of the 60 the selected angular range e.g. O^B^ji, is covered. In 

reference coordinate system 36. The closed contour portion process step 52, the process returns to step 42 and the next 

of the array 34 is defined for each scan line, by an interval scanline is considered. Process steps 42-52 are repeated 

of "x" coordinate values, including a minimum interval until the array is fully scanned, 

endpoint value on the left side of the contour and a FAST COMPUTATION 

maximum interval endpoint value "x r on the right side of the 65 For each angle 0 representing a rotation of the coordinate 
contour. The entire closed contour is thus represented in the system, the global minimum and maximum contour bound- 
reference coordinate system as a sequence of such intervals, ing displacement values are represented by the bounds 
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(5) 



In this sequence of tangents, m e runs through the numerators 

Using simple geometry, the x-displacement of a reference while n e runs through the denominators. Since successive 

system coordinate point (x,y) with respect to axes rotated by elements of these sequences differ by at most one, it requires 

an angle 6 is at most an addition to form a product in equation (4) for a 

given tangent as a function of the corresponding product for 

xe=x cos e+y sin e. (2) the previous tangent. The initial value 2 n is generally con- 
stant for the entire computation. Using this scheme to 

It is preferred to compute equation (2) with integer 15 S enerate 32-sided polygons, the maximum error in angular 

i-p j . , i Q • c—* : approximation occurs at ti/Io and is 0.07 radian. 

operations. To do so, the rotational angle 8 is first appro xi- - .„ A t f , e A A 4 . iL , r 

. , . i * • * „ / „,l jm FIG. 5 illustrates a preferred fast computaUon method for 

mated by an angle having a rational tangent^ m^,, when performing me coordinate system rotation step 48 of FIG. 2. 

the improper expression for the tangent at n/2 is freely * {& ^ iof to im lcmcrjlatiorj of \ hc coordinatc 

written as a shorthand for the corresponding sine/cosine pair. sysl£m rotalioQ step ^ a seq ^ tncc of rational tangents 

Equation (2) can now be written in the form ^ such ^ ^ sequence (5) above) has been selected. In 

a first series of coordinate system rotation steps 54-60, the 

jcq . 1 (ner + mqy). equation (4) components n 6 x J (y) and ngX^y) are computed 

\l ma 2 + ng 2 using the initial n e sequence value 2". For a 32-sided 

„ , polygon formed by 16 rotations of the reference coordinate 

The second step m converting equation (2) to an integer M * ^ 3fi (he yalue rf n fa 2 such ^ (here Me 2 ~ 2=2M6 

computation is based on the observation that computing toe ^ ^ v£jue ^ tfaus 2 „_ 22 _ 4 ^ ^ 

minimum and maximum interval values m equation (1) va]ue ^ Z£rQ m accordance with sequence (5). To calculate 

requires only comparison operations, so that the values in ( } where 2 „ a fa Qm nec to shift lhe value 

equation (3) may be scaled for convenience. Delaying ( } Q Uces ^ cm be efficienfl performed by placing 

multiplication by the (generally irrational) factor x ^ valu6 x (y) for thc rcfer6ncc coordinate 

5 system in one of the processing unit shift registers 20 and 

*e» . . , performing n logical shift left operations. This operation is 

Nm e 2 + *e 2 performed in process step 56 of FIG. 5. Similarly, to 

determine the equation component n e x^y) the reference 

equation (1) becomes 35 system coordinate value x r (y) is placed in a processing unit 

shift register 20 and a logical shift left operation is per- 

, min , , , , (4) formed n times. This operation is performed in step 58 of 

h ~ s * y ww + wqyh FIG 5 afe Wo logical sfaifts left m tne case of a 

mhX 32-sided polygon. The equation (4) value m e y, with the 

r 6 -s Q y {nexjy) + mv), 4Q yahie m ^ Q fe ^ equaJ tQ Z&TQ for ^ fifSt 

rotation angle, x 1 (y,9) is equal to the calculated value 

, _ , , , , nflXi(y), and x_(y,6) is equal to the calculated value n e x_(y). 

where the final products are rounded to the nearest integer. ^ y ^ arc rcturncd ^ ^ 6Q ^ ±c boundi 

An alternative approach is to omit multiplication by the mteryal ya]ues ( e) x ^ Q) - n u 

scale factor s G at this point, leavmg the intervals in propor- 45 If ft ^ dctermillcd in st 54 ±K{ a first rotation has 

uonal umts and performmg inverse transformations on other previously been made> ^ process jumps to step 62 where 

data as required. Although a case can be made for this ^ va]ucs ( Q) and afe dctcrmined 5 sim lc 

alternative m some apphcaUon environments the multiph- ations over a of ^ to r XhuS) 

cation is indicated here in order to simplify later notation. following me tot rotalion using mc sequence value 

Because application of the constant factor s e is delayed in 50 0/2", the next angular rotation employs the m e /n e sequence 

calculating equation (4), exact integer computation of the va i ue w 1 . Because the denominator value remains the 

(proportional) displacement at 0 involves at most three same, the only operation required to map this rotation is to 

multiplies. This allows the implementation of fast compu- a( jd the value m e y to each of the previous x^y.8) and x,(y,6) 

tation methods, such as limiting m/n to direct or inverse values. For m e -l, the value m e y-y. For subsequent 

powers of two, or using multiplication tables when the 55 rotations, the value of m e increases by 1. Thus, x^O) and 

coordinate range is reasonably small. xXy,9) at each subsequent rotation at which ne=2" are 

In order to maintain rotational invariance, at least to the determined by simply adding the value y to the previous 

accuracy afforded by the number of angles chosen, the Xj(y,9) and x,j(y,8) values. The condition m e =0 to 2" is tested 

angles would ideally be spaced equally. But then the integer in process step 62 and the mapped values x^y.S) and x,(y,0) 

values comprising the approximating rational tangents will 60 are returned in step 64 to the process step 44 as bounding 

not have properties that would allow a reduction in the intervals for each new rotation angle. Process steps 62 and 

number of multiphcation operations in equation (4). Instead, 64 are repeated until m 0 =2". At that point, the values x 1 (y,0) 

a set of angles is chosen whose tangents have convenient and x,(y,6) are determined by process steps 66 and 68. 

computational properties, while introducing an acceptably Process step 66 tests to determine that n 0 is in a range 

small error in the even distribution of angles. To divide the 65 between 2" and -2". In this range, the value of n 6 is 

range O=0^ji into 2 W4 " 2 parts, the following sequence of sequentially decremented by 1 until n e =-2". The operation 

me/ne tangents may be used: of decrementing the denominator ^ by 1 changes the values 
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D e x i(y) ^ ^Xy) °y subtracting x 1 and x„ respectively, 
therefrom. The mapped values x 1 (y,8) and x,j(y,9) are 
returned in process step 68. The final sequence of rotations 
is mapped in process steps 70 and 72 by holding the 
denominator value no constant at -2" and decrementing the s 
numerator value m 6 from 0^=2" to m e =l. This sequence is 
achieved by subtracting y from the values Xj^B) and 
x£yfi) f for each rotation. The mapped values are returned in 
process step 72. 

Thus, the mapping operation is easily achieved by an 10 
initial logical shift left operation followed by successive 
addition and subtraction operations. 

A second opportunity for implementation efficiencies 
arises from the fact that there are transition angles associated 
with the interval endpoints as the scan moves between 15 
successive lines. The values of 8 at which an interval 
endpoint minimum or maximum can change are determined 
by the transition angle of that endpoint when proceeding 
from line y to y+1. This function is expressed implicitly, for 
the left end, as 20 

n & ic 1 (v+l)+m e O'+l)<n & r,(y)+m a y, 

and this reduces to 



Similarly for the right end, 

/ne/«e>jf,(y)-xXy+l) 



(6) 25 



(7) 



Process step 48 can be readily modified to test for these 
conditions so that coordinate system rotations that cannot 
yield a Dew interval endpoint minimum or maximum values 
are not tested. FIG. 5 illustrates the importance of transition 
angles in a case where the left interval end displacement 
increases by 1 from line y and y+1 and the right interval end 
displacement decreases by 1 from line y to y+1. Applying 
equation (6), there can be no change in the left endpoint and 
y+1 will not be a new maximum until m e /n e — 11. The right 
end point at y+1 will not be a new maximum until nVne«=l. 
TOPOLOGY 

The term 'topology* is used in a general sense to refer to 
relationships between the above -described convex polygons 
which are set-oriented or metric in nature. Advantageously, 
it is not necessary to compute the actual polygon boundary 
in order to answer questions relating to comparative char- 
acteristics such as containment or proximity. Since the 
interior of the polygon is the intersection of bands, any such 
questions can be expressed equivalently as a conjunction of 
the same questions (or predicate) over the set of component 
bands. 

Abounding polygon "A" of an object may be denoted by 
a set of bands {A^}, where 0 ranges over the set of angles 
that are chosen, and where that set will be the same for all 
objects. Because bands are nothing more than intervals, 
there is an elementary calculus of bands at the same angle, 
equivalent to the calculus of intervals on the line. The 
left-hand and right-hand limits of a band A# may be denoted 
by L(Ae) and R(Aq) respectively. The size of a band is then 

|A e M(AeK(A e ). 

In order to measure the distance between two bands, the 
definition 



30 



35 



A 9 e* Q-maxCLCAflHW, L(B^)~R(A^), 

has the correct properties. If Aq and B e are disjoint, then the 
above value is positive and equal to the distance between 



65 



8 



their closest limits. If the value is negative, then its magni- 
tude measures the extent of overlap. The following are other 
important relations: 

fcy^a^CAa) *x e HR[Aq), 
A^ZBq^L(B q ) (A b ), RtAj*R(b B ), 
A^\B^mzx{L(A e ), L(B Q ))> min(R(A e ), R(B Q )) ], 

The band "axes" thus form a basis for a space of convex 
polygons. The relations listed above can be further extended 
using additional attribute definitions. In the following 
descriptions, roman capital Letters without subscripts denote 
complete polygons over the sent of angles: 

the point can be said to fall within each band of A. For 
inclusion, 

AC0fc*(Ve>4 e V£ e . 

The check is simply that bands of A are contained within the 
corresponding bands of B. Note that for this predicate, and 
for point membership as well, if the result is negative it can 
generally be determined quickly by checking the condition 
in pairs — i.e., a band and its (perhaps approximately) per- 
pendicular band. Two more operations which naturally dis- 
tribute over bands are the following: 

AUB={A e UB e }. 

Based on the measures for bands, one can define 



Width(A) 



6 

' e * 



One can define, as well, the dimension of a polygon along 
any of its axes. Distance generalizes as 



Dist(A,B) = 



AeBBa. 



50 



55 



60 



The maximum over all band distances is the correct measure 
for polygon distance, because no two points on non- 
intersecting convex polygons can be closer, in the Euclidean 
sense, than any band distance. Fortunately, computing abso- 
lute distance is not as important to document analysis as 
testing proximity, which is defined as follows: 

It is not necessary in this case to compute Dist(A,B), but 
only to check the band distance at each angle until either the 
condition is not satisfied or it has been found true for all 
angles. Again, there is generally a rapid exit from the 
procedure for negative results. Thus, a judicious choice of 
the checking order can decrease execution time. 
OBJECT CLASSIFICATION 

A second advantage to defining the bounding polygon in 
terms of fixed directions is that there is then a degree of 
rotational invariance in gauging object structure. It is pos- 
sible to define a line segment, for example, as an object with 
a length -to -width ratio between limits, or with length and 
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width between absolute limits. If only one rectangle is used 
as the bounding figure, only horizontal or vertical lines can 
be detected. When an array of bounding bands are used, it 
is possible, to a tolerance determined by the number of 
angles used, to recognize a line in any orientation. s 

To take a slightly more complex example, consider the 
problem of recognizing dashed lines. Let LMIN be a param- 
eter defining required line thinness, and let MINGAP and 
MAXGAP be parameters delimiting the gap size in a dashed 
line. A recursive step in recognizing dashed lines could 10 
utilize the following comparison; If L is a line segment or 
dashed line, and S is a line segment, then LUS is a dashed 
line if Width(LUS)iLMIN and MINGAP^ Length (LUS)- 
Length (L)-Length (S)^ MAXGAP. 

Recognition of another common graphical shape, the 15 
circle or disk, is specified as condition that all bounding 
bands are, to within some tolerance, of the same size. Other 
approximations to shape may be specified as distributions of 
band size as well. 

These examples illustrate the flexibility inherent in rep- 20 
resenting the bounding polygon as an array of bands in 
accordance with the invention. The array can be seen as a 
basis for any classification which can be done without 
analyzing the interior of the bounding polygon. 
Advantageously, it will be appreciated that weights could be 25 
assigned to the angles chosen, or the angles could be 
parameterized in other ways based for specific applications. 
DRAWING THE POLYGON 

As described above, it is not necessary to calculate the 
vertices of the bounding polygon in order to perform any of 30 
the operations discussed. In graphics applications, however, 
it is often desirable to distinguish an object by drawing its 
bounding polygon. Calculation of the boundary is also 
important in improving estimates of proximity. The prox- 
imity (within A) of two bounding polygons does not dem- 35 
onstrate the point-to-point proximity of the enclosed figures, 
but the error may be reduced considerably by inspecting 
actual picture elements in the interior of a polygon inter- 
section. 

Let A* & be the polygon A extended on all sides by A. 40 
Denote the set of black pels in the interior of a polygon A by 
l b (A). Then a more refined test for the proximity of A and B 
is 



45 



If these conditions hold, then there is point-to-point prox- 
imity along some axis. The accuracy of this test increases 
with the number of polygon sides. 

To calculate the polygon vertices, it is first noted that each 
side of the polygon is a line perpendicular to its 6 axis 50 
through the origin, and at a distance d from the origin. The 
equation of such a line is 



55 



or, in the case under consideration, 

n^-WHey-Vjg (or t^Jq). (6) 

To calculate the vertices, it is only necessary to order the 60 
lines of equation (6), for example, by increasing angle and 
determining the intersections between adjacent lines in the 
(circular) sequence as the polygon vertices. Denoting the 
solution of equation (6) for l e as Lq and for r e as Rq, the 
sequence is 65 

{R Q , . . . , R e< ^>, L^, . . . , Le^, Lq, - - . , L^-cn/^ Rjtfl* 

• • , R 8 <»}. 



The algebra need not be reproduced here, but it will be 
appreciated that the calculations are considerably simplified 
when using the angular sequence (5). 

A class of convex polygons over fixed directions in the 
plane have thus been defined. It has been shown that the 
bounding polygons for all connected components in a binary 
image can be computed efficiently in one pass, and that a 
particular form of the procedure can be implemented using 
only shift and add operations. In addition, a basic algebra 
over the polygons can be used for efficient segmentation, 
low-level discrimination among basic object shapes, and 
recognition of complex objects consisting of disconnected 
components. While preferred embodiments and implemen- 
tations of the invention have been shown and described, it is 
expected that modifications or adaptations thereof will be 
made by persons skilled the art. Accordingly, the scope of 
the invention should not be limited except in accordance 
with the spirit of the following claims and their equivalents. 

I claim: 

1. A computer-implemented pattern recognition method 
for characterizing a closed contour graphics image, includ- 
ing printed characters, using convex-polygon bounding of 
the closed contour image, comprising the steps of: 

selecting an input graphics image containing a closed 

contour character to the recognized: 

inputting a binary array representation of said closed 
contour input graphics image character relative to a 
reference coordinate system; 

generating, in a single pass through said binary array, a 
stored character data output structure representing a 
convex-polygon defining the maximum and mini- 
mum contour bounds of said input graphics image 
character for each of a predetermined number of 
successive, equally sized incremental angular rota- 
tions of said reference coordinate system, wherein 
said contour bounds are defined by parallel lines or 
bands passing through minimum and maximum con- 
tour points on said image relative to each rotation of 
said reference coordinate system; and 

comparing said character data output structure with 
similarly generated character data output structures, 
or with predetermined parameters, in order to rec- 
ognize said input graphics image character. 

2. The method of claim 1 wherein said data output 
structure generating step includes the steps of scanning 
successive lines of said binary array and identifying, for 
each scan line, contour displacement minimum and maxi- 
mum values relative to each rotation of said reference 
coordinate system, and adding said values to said data output 
structure if they represent global minimum and maximum 
values of said input image. 

3. The method of claim 1 wherein said data output 
structure generating step includes the steps of: 

generating and initializing an interval table for represent- 
ing the bounds of the closed contour input image with 
respect to each rotation of said reference coordinate 
system, said interval table defining a set of rotation 
angle approximation values, and a pair of minimum 
and maximum contour displacement values associated 
with each of said rotation angle approximation values; 
filling in said interval table in order to define a contour- 
bounding convex-polygon by scanning successive lines 
of said binary array, and for each scan line: 
identifying successive bounding contour interval pairs 
representing minimum and maximum contour dis- 
placements with respect to each rotation of said 
reference coordinate system; and 
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storing in said interval table one or both of the contour 
minimum and maximum displacement values if they 
exceed the values previously stored at corresponding 
locations in said interval table. 

4. The method of claim 1 wherein said data output 5 
structure generating step includes rotating said reference 
coordinate system by mapping maximum and minimum 
contour displacement values relative to said reference coor- 
dinate system using angular approximation values. 

5. The method of claim 4 wherein said coordinate system 
rotating step includes determining the maximum and mini- 10 
mum contour displacement values "x a ,y" and "x„y" relative 

to said reference coordinate system directly from said binary 
array, and generating mapped minimum and maximum 
contour displacement values 1 B " and "r e " by modifying said 
"x 3 ,y" and "x„y" coordinates in accordance with said angu- 15 
lar approximation values. 

6. The method of claim 5 wherein said angular approxi- 
mation values are a sequence of rational tangents of a form 
me/rie, wherein successive m e and successive n 0 values 
differ by at most one. 

7. The method of claim 6 wherein said rational tangents 20 
m/n are selected to divide the angular range 0<8<n into 2 n+2 
parts in accordance with the predetermined sequence: 

{0/2", . . . , 272", 272"-l, . . . , 270, 27-1, . . - , 27-2", 
2"-l/-2 w , . . . , 1/-2"}. 

8. The method of claim 7 wherein said coordinate system 25 
rotating step includes determining the mapped minimum and 
maximum contour displacement values "l e " and "r e " in 
accordance with the relationships l e -n 0 x J (y)+m e y, and 
r e *n ebeJ (y)+m e y, over the predetermined sequence of m/n 
rational tangent values. 

9. The method of claim 8 wherein said coordinate system 30 
rotating step includes, for each scan line "y" of the binary 
array in said reference coordinate system, the step of map- 
ping the minimum and maximum contour displacement 
values "Xj,y" and "x r ,y M to the rotated coordinate system 
values "l e " and "r e " by placing the values "x r " and "y" 35 
in data storage register locations and performing shift and 
add operations on said values. 

10. The method of claim 9 wherein said coordinate system 
rotating step includes a first step of placing the and "x r " 
displacement values in first and second data storage register 40 
locations, and placing the scan line value "y" in a third data 
storage register location, a second step of determining the 
values "viqX^ and "n e x r " for the start of the predetermined 
m/n tangent sequence by performing a logical shift left on 
the "x/* and x/' values, and a third step of determining the 45 
mapped interval pair values "^=(^1^+0^)" and "r e = 
(neXj/t-mey)" by successively adding and subtracting the 
values "y", "x/* and "x r " in accordance with the change in 

m and n over the predetermined sequence of m/n tangent 
values. 

11. The method of claim 10 further including the step of 50 
limiting the sequence of tangent m/n values over which said 
reference coordinate system is rotated upon determination 
that the difference between the and "x" values for 
successive scan lines of the binary array are greater than or 
less than, respectively, the m/n tangent value. 55 

12. In a programmable data processing apparatus includ- 
ing a central processing unit, a data input system, a data 
output system, and one or more data storage resources, a 
pattern recognition system employing character bounding 
for evaluating closed contour images, comprising: 60 

means for inputting a binary array representation of one or 
more closed contour objects to be evaluated with 
respect to a reference coordinate system, said closed 
contour objects representing an input graphics image 
containing closed contour characters to be recognized; 55 

means for generating from said binary array representa- 
tion one or more stored convex-polygon representa- 



tions of said closed contour objects of a form that 
includes a set of polygon bands oriented relative to said 
reference coordinate system and to a predetermined 
number of successive, equally sized incremental angu- 
lar rotations of said reference coordinate system, said 
bands representing a convex-polygon defining the 
maximum and minimum contour bounds of said closed 
contour objects with respect to each of said angular 
rotations of said reference coordinate systems, wherein 
said polygon bands include pairs of minimum and 
maximum contour displacement values representing 
the bounds of the closed contour objects for each 
rotation of said reference coordinate system; and 
means for comparing said convex-polygon representa- 
tions with similarly generated convex-polygon 
representations, or with predetermined parameters, in 
order to recognize said closed contour objects. 

13. The apparatus of claim 12 further including means for 
comparing plural convex-polygon representations of said 
closed contour objects in order to isolate said objects and 
approximations of their shape. 

14. The apparatus of claim 13 wherein said means for 
generating convex-polygon representations include a first 
means for scanning said binary array representation, a 
second means for determining the contour bands relative to 
said reference coordinate system, and a third means for 
incrementally rotating said reference system contour bands 
and adjusting them to bound the closed contour objects in 
the rotated directions. 

15. The apparatus of claim 14 wherein said third means 
includes fourth means for mapping minimum and maximum 
contour displacement values from the reference coordinate 
system through successive rotations thereof by successively 
incrementing and decrementing the reference coordinate 
system contour displacement values in accordance with a 
predetermined sequence. 

16. The apparatus of claim 15 wherein said fourth means 
includes fifth means for limiting said predetermined 
sequence upon determination that the change in minimum 
and maximum contour displacement values between suc- 
cessive scan lines of the binary array has reached a critical 
threshold. 

17. The apparatus of claim 16 wherein said critical 
threshold is the approximate angle by which the reference 
coordinate system has been rotated. 

18. A computer-implemented pattern recognition method 
for characterizing a closed contour graphics object using 
convex-polygon bounding of the closed contour object, 
comprising the steps of: 

inputting a binary array representation of said closed 
contour object relative to a reference coordinate sys- 
tem; 

generating, in a single pass through said binary array, a 
stored data output structure representing a convex- 
polygon defining the maximum and minimum contour 
bounds of said input closed contour object relative to 
successive, uniform rotations of said reference coordi- 
nate system, said contour bounds being determined by 
mapping pairs of interval endpoint values of the refer- 
ence coordinate system by successively incrementing 
and decrementing the "x" and "y" coordinate values 
thereof in accordance with a predetermined sequence; 
and 

comparing said data output structure with similarly gen- 
erated data output structures, or with predetermined 
parameters, in order to recognize said closed contour 
object. 
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